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Abstract 

In this paper we consider the soft processes at LHC energies in the framework of the Con- 
stituent Quark Model (CQM). We show that this rather naive model is able to describe all 
available soft data at lower energies and to predict the behavior of the total cross section, elas- 
tic and diffractive cross sections at the LHC energy. It turns out that the "input" Pomeron, 
which has been used in this approach, has parameters that are close to so called "hard" Pomeron 
with rather large intercept A ~ 0.12 and small value of the slope a'p « 0.08Gey~^. We show 
that the elastic amplitude has a minimum at impact parameter 6 = and a maximum at 
b ~ 2GeV^'^. Such a behavior is a result of overlapping the parton clouds that belong to 
different quarks in the hadron. 
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1 Introduction 



One of the most challenging problems of QCD is to find the correct degrees of freedom for high 
energy " soft" interactions. The question is what set of quantum numbers diagonalizes the interaction 
matrix at high energies. The Constituent Quark Model (CQM) pp is one of the models which can 
be a good candidate for a correct descriptions of the "soft" interactions, see also j2]. In this model 
the constituent quarks play roles of the correct degrees of freedom for high energy QCD and the 
structure of a hadron is characterized by two radii: the proper size of the constituent quark {Rq) 
and the typical distance between two constituent quarks in a hadron {R). The main assumption is 
that R :$> Rq. 

In spite of the fact that this model looks rather naive, it is supported by the two sets of the 
experimental data, namely, CDF double parton cross section at the Tevatron, [Sj, and HERA data 
on inclusive diffraction production with nucleon excitation, jHEI- In our paper , [HI, we examined 
these data and found, that the CQM model describes a lot of "soft" data in the first approximation, 
see also [Tj. The radius of the constituent quark, which was found in Ref. ^6], turned out to be small: 
Rquark ~ O-^ ^ 0.2GeV~^. However, this radius depends on energy (at least logarithmically as 
R = R^uark + 4ap \n{W /Wq)), and the possible scenario is that at the LHC energy this radius 
becomes compatible with the distance between constituent quarks {R) (see Fig. (^). Therefore we 
could expect a new physics at the LHC in such an approach. 



W < LHU energy W > LHU energy 




Figure 1: The protoTi structure at different energies in CQM. 

In this paper we are going to develop a systematic approach to high energy scattering in CQM. 
To our surprise we found that only the simplest diagrams of this model has been discussed in details, 
namely the diagrams which include only interaction between a pair of quarks. In our model we include 
all possible quark interactions, considering an eikonal approximation for the scattering amplitude of 
two colliding quarks, (see Fig. (j21)). In this case, the first contribution to the elastic amplitude will 
be simply nine interactions between constituent quarks in the protons. Further contributions to the 
amplitude will include the other interaction between quarks (see Fig. 0). Taking into account all 
possible configurations of quarks we will obtain amplitude which contains nine different contributions 
with the alternating signs. Such a structure of the answer is similar to the scattering amplitude of light 
nuclei (tritium-tritium scattering). We will show, that the result satisfies the unitarity constraint, 
when we consider interactions between different configurations of the constituent quarks. The effects 
of interaction of several pairs of quarks is especially important at high energies. If we will consider 
our amplitude at asymptotically very high energy, where we will may replace eikonalized amplitude 
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by a step function, the different parts of the amphtude will cancel each other and only the last 
diagram, in which all quarks of the projectile interact with all quarks of the target, will survive. This 
last diagram will give unitarized Froussart-like answer. In our estimates it turns out, that already 
at energies a/s = 1855 GeV we need to consider the interaction of the five pairs of quarks in the 
protons. At the LHC energies, the interaction of the seven quark pairs is essential. This structure 
of the interaction changes not only the high energy dependence of the scattering amplitude but also 
leads to a quite different impact parameter dependence of the answer. In the case of these multi- 
Pomeron exchanges between the different pairs of quarks, the amplitude has a minimum at low b. 
Such change is very important since it could affect the behavior of the slope both as function of 
energy and as function of momentum transfer. 

Another interesting problem, addressed in this paper, is the values of parameters of the "initial" 
Pomeron. Indeed, it is widely believed, that the " soft" Pomeron is originated by the non-perturbative 
QCD contributions, which are out of theoretical control at the moment. Everything that we know 
about 'soft" Pomeron is a mixture of our phenomenological knowledge with the general theoretical 
ideas on the properties of non-perturbative QCD contributions (see Refs. jHl El UHl HH 1121 UHl UH 
IT5| IT6| IT7\ ITSj IT9] ) . The question is the following, can we obtain the well known phenomenological 
Pomeron using as "initial" the Pomeron with the large intercept and small slope, i.e so called "hard" 
Pomeron, jSOj ? Does the "hard" Pomeron play any role in the "soft" interactions? There exist two 
different points of view on this question. The first one is that we need to introduce separately two 
objects, "soft" and "hard" Pomerons, which have different properties and contribute differently, each 
in a different kinematic region, see for example jHl El 1^ ^| and references therein. Another point of 
view, see for example j2] , is that non-perturbative physics comes in our calculations only in the form 
of the boundary and/or initial conditions and the "soft" Pomeron arises as a result of unitarization 
effects and self- interactions of the "hard" Pomerons in the amplitude. Our model may help to clarify 
the situation. Taking into account all effects of the unitarization, eikonalization and accounting all 
interactions between the different configurations of quarks, we will fit the experimental data. The fit 
will determine the parameters of the "initial" Pomeron , such as its intercept and slope, as well as 
the radius of the constituent quark. We show in this paper, that the parameters of the initial "soft" 
Pomeron are close to the parameters of the " hard" Pomeron and quite different from the parameters 
of the Donnachie-Landshoff Pomeron [22} , which is usually considered as a typical " soft" Pomeron. 

The structure of the paper is as follows. In Section 2 we discuss in more details our approach 
and methods of the calculations. In Section 3 we apply our model to the p-p data and fit the 
experimental data in order to find numerical values of the parameters of the "initial" Pomeron. 
Section 4 is dedicated to the elastic amplitude as a function of the impact parameter. In this section 
we also consider the different contributions to the elastic amplitude due to interactions of different 
numbers of quark-quark pairs. In Section 5 we calculate the survival probability (SP) of the exclusive 
hard processes in p-p scattering . The cross section of the diffractive dissociation process is calculated 
in the Section 6. The last section. Conclusion, contains the main results of the paper as well as the 
discussion of a future work in the proposed direction. 

2 Proton-proton scattering in the Pomeron approach 

The key ingredient of the CQM is the quark-quark scattering amplitude. Considering this model, 
we need to determine the form of the single Pomeron exchange between two quarks. The next step 
will be the eikonalization of the single scattering amplitude, that means a replacement of the single 
Pomeron exchange in the scattering of the particular pair of quarks by the eikonal amplitude. The 
third step in our calculation will be the consideration of the interactions between all possible quark 
configurations in colliding protons, see Fig. El For this last step we need to know the wave function 
of the quarks in the proton and the vertices of the Pomerons-quarks interactions. Only after the 
determination of the wave function and vertices we will be able to calculate the diagrams for the 
quark-quark interactions. 
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Now let us consider these problems step by step. 



2.1 Quark-quark interactions 

We determine the amplitude for q-q and p-p scattering in the impact parameter representation. 
In this case the "soft" Pomeron exchange for the interaction of the pair of the quarks (Pomeron 
propagator) has the following form (see more about soft Pomeron in Ref.|18 [ l22 | I23j: 

n,-,iY,b) = ^oe^"-'-^, (1) 

here, F = In GeV'^) is the rapidity of the process, b is the impact parameter of the process, A 
is the intercept of the 'initial', input Pomeron and 

R"^ = 8Rl + Aa'p ln(s/so) . (2) 

Here Rq is the squared radius of the constituent quark and a'p is the slope of the input Pomeron 

trajectory. The numerical values of the cxo , A , R^ and a'p we will find fitting data for p-p scattering. 
The eikonal amplitude, which is a "main" ingredient in our calculations, we determine as follows: 

Pg-g{Y,b) = 1 - e-^'-'(^'^)/2 ^ 

see Fig. El where Pq-q is the imaginary part of the quark-quark scattering amplitude at high energy. 
Below , discussing the single quark-quark interaction amplitude, we will mean only the amplitude 
given by Eq. Q. The Pomeron of Eq. we will consider only as the "input", initial Pomeron of 
our problem. 




Figure 2: The Pomeron exchange and eikonalized quark-quark amplitudes. 



2.2 The model of the proton 

In order to take into account all possible configurations of the interacting pairs of the quarks, see 
Fig. 121 and all figures in the Appendices, we need to know the analytical expressions for the vertices 
of the quark-Pomeron interactions. It is clear, that we need to calculate only three types of such 
vertices, see Fig. EJ-Fig. where there are one, two or three groups of the Pomerons attached to the 
one, two or three quarks. In order to calculate these vertices we need to know the wave functions of 
the constituent quarks inside a proton. We use a very simple Gaussian model for this wave function, 
which corresponds to the oscillatory potential between pair of quarks in a proton. In this model we 
have ( see Pl]): 
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Figure 3: The three first orders of the p-p elastic amplitude in the CQM. 



where the constant a is related to the electromagnetic radius of the proton: 

a = l/i?Lcir ^ OMGeV . 
For the diagram of Fig. EJ we have: 

r 2 

Viiq) = / |^(Xi,X2,X3)p5(Xi+f2 + ^3)e*«^^^ = 6"^ 



(4) 



(5) 



(6) 



= e 



-qV6a 




Figure 4: The one Pomeron vertex in quark-quark interaction. 

We determine the impact factor for two groups of Pomerons attached to two different quarks (see 
Fig. ED : 
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V2{qi,q2) 



2 2 ^ ^ 



(7) 



2 2 




= P 6a 6a 6a 



2_ + 



Figure 5: Two Pomeron vertex in quark-quark interaction. 
The last impact factor is the vertex which is shown in Fig. IHl We have for this vertex : 



(?2-g3)^ (<ri-?2/2-g3/2)^ 

g 8 Q 6 Q 



In all three cases the total transferred momentum of the diagrams is defined as the sum of transverse 
momenta of all Pomerons (qi): k = J2 qi ■ 




-(^2-^3)' - (qi-q2/2-q3/2)' 



= e 



8a 



6a 



Figure 6: The three Pomeron vertex in q-q interaction. 



2.3 The elastic amplitude of p-p scattering 

We have all ingredients for the calculation of the elastic amplitude. In the CQM we can write the 
amplitude as the sum of the amplitudes with the different number of interacting quark pairs: 

A{S, h) = Aipair{s, h) - A2pairs{s, h) + A^pairsis, 6) + ... + A^pairs{s, h). (9) 

Here the amplitude for one quark pair interaction Aipair is equal to Pg-q which is defined by Eq. dHj). 
The maximum number of possible quark pair interactions in the amplitude is 9. The three first 
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orders of the possible configurations of interactions are sliown in Fig. IHl and the calculations of all 
other orders are presented in the Appendix B. The amplitude defined in a such way incorporates 
unitarity by construction. Indeed, checking the expressions for the different term of the amplitude, 
that are written in Appendix B, it is easy to see, that at asymptotically high energies only the last 
term of Eq. will survive, giving the Froussart-like answer for the whole amplitude. 

To calculate the contribution of the first diagram of Fig. El i.e. the first term of the r.h.s. of 
Eq. 0. we make Fourier transform of the vertex of Eq. (jHl) from momentum to impact parameter 
space: 



47r2 27r 

Using this vertex we obtain the contribution to the elastic amplitude from the one Pomeron exchange 
Aig_q{s, b) , which is : 

A,,aUs, b) = 9^^J dH, J dH, e-^(^^-^^)^ P,_,(y, b,) e-^ , (11) 

where the first coefficient in Eq. (jllj) is the total number of q-q interactions in this order. Using more 
complicated vertices, which are given by Eq. ((7|) and Eq. (jSJ, we calculate all terms that contribute to 
the elastic amplitude, i.e. terms on the r.h.s. of Eq. Q. In Appendices A,B the resulting expressions 
for the amplitude A{s, b) are written as well as the examples of the calculations of the diagrams 
with the different numbers of interacting pairs of quarks. 



2.4 The total, elastic cross sections and the elastic slope Bei of the 
proton-proton interaction 

The expression for the amplitude (Appendix B), is determined in the impact parameter space and 
now we can easily calculate the different cross sections for the processes of interest. 
(Jtot- The total cross section in impact parameter representation is simply 

atotis) = 2 I dHlmA{s,b). (12) 

Fitting the experimental data at low energies we also add to this cross section the contribution of 
the secondary Reggeons, aj^l^ , see j2S] • 

(Tel. The elastic cross section in the same framework is equal: 

a,i{s) = J dH\A{s,b)\^ = J d^b{ImA{s,b)f. (13) 

where we assumed that the quark-quark amplitude is mostly imaginary as it follows from the eikonal 
approximation. 

Bel. We consider only the first term of the elastic amplitude, i.e Aiq_q{s,b) . In this case we 
have for the slope B^i : 



9a 
4 



^JdH,JdH2jdH b' e-^(^-^i)' P,_,(y, 62) e- 



3n 
2 



(fel-fe2)2 



2/ dHPq^q{Y,b) 



1 crlr ^ I dH b^ Pq_q{Y,b) 



a o- 



Iq-q 
tot 



\q~q 



(14) 
(15) 



tot 



Here 



Iq-q 
tot 



2 J dHPq,q{Y,b). 
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We add the contribution of the secondary Reggeons to this expression , which is needed to be taken 
into account at low energies: 



^el = „ Ig-g , Reg Klectr + „ Ig-g , Reg + 2 (a^ " a^p) ^^[s/So) j^^g " (16) 

The third term of the r.h.s of Eq. (fT^ contains the contribution of the secondary Reggeons, which we 
do not consider in our model. The secondary Reggeons we account only on the level of proton-proton 
scattering. Therefore, the parameters which we take for this contribution have a pure phenomeno- 
logical origin. We take for the slope of secondary Reggeons a^i = 0.86 GeV'"^ and for the slope of 
phenomenological "soft" Pomeron ^ a^p = 0.25 GeV~^ at a/sq = 9GeV in the r.h.s of Eq. (fTTIjl . 
see |22l EHl ES] • Generalizing this expression to the case of full elastic amphtude ( see Eq. ©) we 
obtain at low energies: 

P 9 air p2 , 9/ dHb^ P,^,iY,b) 

- ■ R^ ^electr + \ RTg + 1^0 

(^tot + CTtot a tot + (^tot 

,,,,,, a^^^ J dHb^ ImAi(b,s) , , 

+ 2 (a^ - asp) In {sjs^) _ + E ' , Rig > (18) 



(^tot + 0"tot i=2 '^tot + Cto^ 



whereas at high energy we have 



R 9^p2 , 9/ ^^&fe^P,_,(y,&) , ^ \ dHb^ ImA,{b,s) 

^el = ^electr H ^ • lJ-9) 

It is also interesting to calculate the elastic cross section using simple expression for the elastic 
cross section which is obtained in the model with the phenomenological "soft" Pomeron. Indeed, in 
this case we calculate the elastic cross section using the following formula: 

Ms) = T^. (20) 
16 vr Bel 

This expression is correct only in the case when one "soft" Pomeron is considered. It will be inter- 
esting to compare the calculations of elastic cross section given by Eq. (|T!^ with the calculations of 
Eq. pUjl . Indeed, in this case we check the possibility to reproduce the simple result of Eq. by the 
theory where many eikonalized Pomeron exchanges are taken into account. So, in the next section 
we will perform the data fitting and will make the calculations using both expressions, Eq. (fT^ and 
Eq. ()20|) in order to show the importance of many quark pairs interaction. 



3 The proton-proton scattering data and the parameters of 
the input Pomeron 

Now, we are able to apply our model to the p-p interactions and, fitting the experimental observables, 
we will extract the values of the parameters of our input Pomeron ( see Eq. (^). In Fig.EJwe present 
the plots for the total cross section, elastic cross section and elastic slope in W = 23 — 1855 GeV 
energy range. We perform all calculations numerically ^ using the formulae of the previous subsection 
with the amplitude written in the Appendix B. There are two different plots which we present for 
elastic cross section using definitions Eq. ()13p and Eq. ((201) • The solid line represents the elastic cross 
section given by Eq. ()2n|l and dashed line represents calculations performed with Eq. ()13|1 . 

^The slope a^p is the slope of the phenomenological Pomeron in proton-proton scattering , see [321123 EH] j and 
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Figure 7: The plots for the total cross section, Fig. elastic cross section , Fig. (0j-c, and elastic 

slope, Fig. ^-h. The experimental data are from 125] and references therein. 
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Figure 8: The plots for the total cross section, Fig. ^-a, elastic cross section , Fig. (0)-c, and elastic 
slope, Fig. ^-h, at high energies. 
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From these plots we see, that the model describes the experimental data quite well. The param- 
eters of the single Pomeron of Eq. (0) extracted from the data fitting are the following: 



• the slope of the input Pomeron trajectory: 



a = 0.08 Gey 



-2 



(21) 



• input Pomeron's intercept: 



A 



0.118: 



(22) 



• the value of cross section for the input Pomeron at sq = 1 GeV^ : 

ao = 6.3 GeV'^; 



(23) 



• the radius of the constituent quark: 



'quark 



0.16 GeV 



-2 



(24) 



With these parameters we extrapolate our calculations for the cross sections and slope at higher 
energies. The resulting plots are shown in Fig. |H1 

It should be stressed that the above parameters are quite different from the Donnachie-Landshoff 
Pomeron [22] • The Pomeron intercept, which we obtained, is higher and the slope is much lower 
than the intercept and slope of the D-L Pomeron. We may conclude, that such small value of the 
Pomeron slope indicates that our input Pomeron can have a "hard" origin. 

4 Behavior of the elastic amphtude of p-p scattering in our 
model 

The elastic amplitude in our model has a sufficiently complex structure. It contains contributions 
from different configurations for the interactions of quarks inside the protons. We found all possible 
terms which contribute to the elastic amplitude, see Appendix B, and, therefore, we can find out 
which terms in elastic amplitude are important at different energies. It turns out, that even at low 
energies the two quark pairs interactions, i.e. interactions between two quarks in one proton with two 
quarks in another proton, are not negligible. At Tevatron energy, ^/s = 1855 GeV , we already need 
to take into account the contributions from the five pairs of interacting quarks. At energy of order 
of the LHC energy, ^/s = 15000 GeV , the contribution of the interaction of the seven quark pairs 
is valuable. At this energy the contribution of two quark pair interaction exceeds the contribution 
of one quark pair interaction, see Fig. IHl This is a signal that the parton clouds start to overlap, 
leading to the picture of Fig. ((1} . 

From Fig. Elwe also see, that even at ^/s = 1855 Gel^ energy, the one quark pair contribution to 
the elastic amplitude exceeds unity . The unitarization of the amplitude in this case is achieved not 
by eikonalization of the quark-quark interaction but by including more complicated configurations of 
the quarks inside of the protons in the elastic amplitude. Therefore, the form of impact parameter 
dependence of the amplitude turns out to be different from the usual Gaussian one. Indeed, the 
contribution of two quark pair interactions, which have a negative sign in the amplitude is equal 
or larger than contribution of one quark pair interaction. At the same time the one quark pair 
interaction is wider in the impact parameter space, see again Fig. IHl The contributions of all terms 
in the elastic amplitude lead, therefore, to the situation where the maximum of the elastic amplitude 
moves from the zero impact parameter to the impact parameter b ^ 2 GeV~^ , see Fig. ^1 This 
effect reflects very simple physics. At high energy at small impact parameter the multi quark pair 
interactions are important and elastic production becomes mostly peripheral. 

it is not related to the slope ap of our "initial" Pomeron in the quark-quark scattering. 
^The Fortran code can be sent after the request (email:sergb(8)mail. desy.de) 
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Figure 9: The contributions to the elastic amplitude from the interactions of different numbers of 
quark pairs at different energies. 
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Figure 10: The elastic amplitude of the p-p scattering as function of impact parameter 
at different energies. 
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5 Survival probability (SP) of the 'hard' processes in p-p 
scattering 

In this section we consider the calculation of the survival probability for the process:p + p — *■ p + 
[LRG] + dijet + [LRG] + p, where p and p are the colliding protons and LRG is the large rapidity 
gap. In this process there are two large rapidity gaps - the intervals of rapidity without secondary 
hadrons. The cross sections of such processes are small in comparison with the inclusive production or 
in other word, in comparison with the process where no rapidity gap is selected. The ratio of these 
two processes, exclusive and inclusive (or the process without LRG), we call survival probability 
of large rapidity gap (see Ref. IZBij). In the simple case of the eikonal approach to proton-proton 
interaction, the survival probability is determined by a simple formula (see R,ef . |26[ l?7] ) : 

^2 ^ / d?hA{s,h) ahard{b) ^^5) 
J d'^b ahard{h) 

where 

A{s,h) = e-""^'''^ (26) 

is the probability that no inelastic interaction between the scattered hadrons has happened at energy 
^/s and impact parameter b . 

Using a simple Gaussian parameterization for ahardip) , namely, 

_ fc^ 
e ^ 

O'hardib) = CThard ^_ ' i'^'^) 

Rjj = 8 Gel/-' (2^ 



32 



with 

we reduce Eq. fE3|) to the form 

= [ dHA{s,b)e~^ . (29) 

vr Rfj J 

The difference of the calculation of the SP in our model from the calculations above is that 
we consider as principle degrees of freedom the constituent quarks but not protons. Therefore, we 
need to determine the expression for SP in terms of the interacting quark pairs. We begin with the 
discussion of the cross section for the exclusive "hard" production in the case of interaction of only 
one pair of quarks. For this process we have: 

O'hardib) = Ohard ( ^lpair(s, &)) , (30) 

with the amplitude Aipair{s,b), calculated from Eq. (|TT|l with he following replacement 



2R 



P,-,(r, 62) - ^-—^P,-,{Y, 62) , (31) 



where 



P,-,{YM) = 1 - P,-,iY,b2) = e-^^-^(^'^^)/2 (32) 

in Eq. (jHTj) means that t no inelastic interactions occur between the considered pair of quarks. This 
pair of quarks, which produces "hard" jet, interacts elastically. The amplitude of this process is 
illustrated in Fig.[TTJ In expression for the A{s, b), we need to introduce a new "hard" radius Rq^h ; 
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Figure 11: The first order amplitude of the LEG "hard" process. 

which related to the "hard" process on the quark level. The numerical value of Rq^h 

the following way. For simplest process of "hard" production without any "soft" rescattering, the 

answer will be the same for any model : 

e"^ = f dH, [ dH.e-'^^'-''^' e"^ (33) 

Atc-^ J J 

Indeed, we can see, that using this equality in Eq. (jHUjl . we reproduce the expression given by Eq. (j2Zj). 
From the Eq. (jSSl), we obtain the value of Rq^h '■ 

Rl_H ^ 1.1 GeV-^ . (34) 

Generalizing the procedure given by Eq. ()3H1 we obtain the answer for SP amplitude A{s, b) in the all 
orders. Indeed, any term of the amplitude given in Appendix B, see Eq. ()B.1|) . have the integration 
over the product of the eikonalized amplitudes of the interacting quark pairs Pq^qipi) ...Pq_q{hj) . 
When we calculate the SP, we should take into account the possibility that the "hard" process can 
occur in the interaction of any pair of quarks. Therefore, for the pair with the "hard" process 
the replacement of Eq. (p?T|) should be performed, whereas other pairs of quarks will still interact 
elastically. It means, that for any integral in expansion of Eq. (jB.lj) of elastic amplitude, we will 
make the following replacement for any term of Pq^qipi) ... Pq_q(hj) - type : 



Pq.q{h) ... Pq.q{h^) ^ \ ^ P,~M ... P,_,(6,) (35) 

for each Pq^qipk) in the chain Pq^qipi) ... Pq^qihj) . This procedure is illustrated by Fig. [T21for the 
case of the interaction of two pairs of quarks. 

Performing these replacements and summing all terms, we obtain a new amplitude A(s, h) ( see 
Eq. ()B.1|) ). In Appendix C we show the example of this procedure for the A^pairs{,s,h) term of the 

elastic amplitude A{s,h) . With this amplitude, A(s, 6) , we determine the "hard" cross section as 
follows: 

(Thard (X. j d^h (i(s,6)) , (36) 

where the first term of this expansion has the form: 



-lard OC J dH ('^J dH, J dH, e-^(-^^)^ Pq^qiY, h,) e-'-f (^V^^)^ ] (37) 
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Figure 12: The general structure of the amplitude of the LEG "hard" process for the 
case of interaction of two quark pairs. 



SP{^/s = 23GeV) 


5P ( = 1855 GeV 


SP{./s = 15000 GeV) 


SP{^/s = 18000 GeV) 


0.23 


0.052 


0.0175 


0.0155 



Table 1: The survival probability factor as a function of i/i 



Finally, the answer for the survival probability factor in all orders of the expansion of Eq. (jB.lj) looks 
as follows: 



J d'b {A{s,b 



2 ■ 



(3^ 



For example, there is the first term, which contributes to the S'^ : 



JdH dH,f dH.e-'^^'-'^^' e A-,(F,62)e-^(^- 



■b2? 



Si 



JdH l^JdH^JdH,e 



(39) 



The values of the survival probability factor, S"^ , calculated in this approach, are shown in the 
Table H We see, that these values are close to the values of survival probability for the central 
diffraction process calculated in Ref. j2H|- We also performed the calculation of the integrant in 
the numerator of Eq. (PH|) , which is proportional to the amplitude of the "hard" central diffraction 
process. The plots of this integrant as a function of impact parameter at different energies are 
presented in Fig. El Considering these graphs we conclude, that the main contribution to the 
amplitude of the central diffraction "hard" production comes from the non-central values of impact 
parameter and it is almost zero for the central impact parameter at high energies. 
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0.06 
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b(GeV'^) 
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b(GeV'^) 



Fig. (O-a 



Fig. (P|)-b 



Figure 13: The impact parameter behavior of the integrant in the numerator of Eq. at energies 
W = 1855 GeV and W = 18000 GeV . 



6 The cross section of the diffraction dissociation 

The elastic amplitude of the considered model accounts many Pomeron exchanges between different 
pairs of quarks. Now, we want to calculate the cross section of the diffraction dissociation processes 
taking into account all these processes. The diffractive dissociation (DD) processes for us are all 
processes where no particle have been produced in the rapidity region between two scattered protons. 
From the unitarity constraint we have: 

O'tot = O-el + (Jinel + (TdD ■ (40) 

Therefore: 

(^DD = CTtot - CTel " CTinel- (41) 

From our previous calculations we know the values of aei and cxtot ■ So, in order to calculate o"£)£) 
we need to calculate the value of amei ■ The calculation of amei is pretty simple. In the first order, 
where 

18 J dHP,_,{Y,b), (42) 



9 / dHP;j^{Y,b), (43) 



Pq-^giY,b) = 1 - e-^^'-'(^'^^ (44) 



a 



total 



2 / d'^bAipairis,b) 



for cr,„p/ we have: 



with 



a,nel = J dHA\lt{s,b) 
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The calculation of amei for all possible interactions between the pairs of the quarks in the protons 
we perform using the following simple receipt. In each expression for atot 

cTtotai = 2 I dHA{s,b), (45) 

we make replacement: 

P,.,{Y, 6.)...P,„,(F, b,) i P;-J(F, h)...Pr4{Y, 6,) , (46) 

with the given by Eq. ()44j) . After the calculations of (Xmei with the help of Eq. ()4H|1 we obtain 

the value of a or, . The result of calculations for the energy range W = 100 — 3000 GeV is shown 
in Fig. 1131 as well as the sum of elastic and inelastic cross sections in comparison with the total cross 
section. 




W(GeV) W(GeV) 
Fig. (Hl-a Fig. (Pj) -b 

Figure 14: The sum of elastic and inelastic cross sections in comparison with the total cross section 
and DD cross section as the function ofW = i/i. 

For the energies higher than W = 3000 GeV , the cross section of diffractive dissociation process 
is almost zero, see Table |2l as it must be at high energies when we approach the "black disc" 
limit. Nevertheless, it is important to notice, that in calculations of DD processes, we neglected the 
diffraction dissociation processes with the large mass production. Such processes may be described 
if we introduce in the model the triple Pomeron vertex, see We did not consider this vertex in 
our calculations, therefore, our result for the cross section of the DD processes, which is the sum of 
single diffraction and double diffraction processes, is smaller than it could be in the case when the 
triple Pomeron vertex is included. 
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o^ddW^ — 7(J(J(JGev/ J 


CJ^DDlv^ = lUUUU GeK j 


c^ddIv-^ = 15UUUGel/j 


<^ddW-^ — loUUUGel/j 


z.o mb 


U.4 mb 


0.35 mb 


U.o mb 



Table 2: T/ie diffractive dissociation cross section at high energies as a function of y/s. 



7 Conclusion 

In this paper we consider the proton-proton interactions in the framework of the Constituent Quark 
Model (CQM) at high energies. The typical "soft" process of p-p scattering we described taking 
into account all interactions between pairs of the quarks in the protons, modeling the quark-quark 
interaction by eikonal formula. It turns out , that in this model the interactions of one, two or 
three quarks from one proton with two or three quarks from the second proton surprisingly become 
essential even at low energies. Indeed, if we look at Fig. IHl we see, that even at energy y/s = 23 GeV 
the contribution from the interaction of two pairs of quarks is approximately one third from the 
contribution of one pair. Of course, at higher energies the contributions of more interacting quark 
pairs became to be more important. For example, at LHC energy the interactions of one and two 
pairs of quarks equally contribute to the elastic amplitude and at this energy we need to account the 
contribution of seven quark pairs. This interaction picture leads to the very natural scenario of the 
unitarization. At high energies the contributions from the interaction of one and two pairs of quarks 
cancel each other, and final amplitude at given impact parameter is smaller than one, see Fig. 
We see, that in this case the unitarization is achieved by exploring internal structure of the proton 
rather than details of interactions. However, we should stress, that we assume that the quark-quark 
interactions is described by the eikonal approach. This specific mechanism manifests itself in the 
form of the b-dependence of the amplitude, which is quite different from the usual b dependence. 
Our amplitude has maximum at 6 ~ 2 GeV~^ at high energy, that means that elastic production 
is mostly peripheral at high energies. 

Another result , observed in the present model and related to the unitarization of the amplitude, 
is the relatively small value of the amplitude at zero impact parameter. Indeed, even at energy 
y/s = 18000 GeV^ at 6 = the amplitude is not close to one. The interactions are still "grey" and 
not "black". Many Pomeron interactions between different quark configurations in both protons lead 
to the "grey " picture in spite of the large contributions of the one, two or three interacting quark 
pairs. 

The considered model fits the experimental data pretty well, see Fig. [3 Using the parameters of 
the "input" Pomeron , which we obtained through the data fitting, we can predict the values of the 
cross sections at high energies. Doing so, we did not find the deviation of the behavior of the elastic 
slope, Bf^i , from the simple linear parameterization, see plot of Fig. |H| The explanation of this effect 
is the following. In spite of the hope, that the constituent quarks will have strong overlapping at 
high energies and this effect will change the energy behavior of the Bei , it actually does not happen. 
The value of the size of constituent quark obtained from the data fit is small, Rg^ark ~ 0-16 GeV'^ . 
The slope of the initial quark-quark amplitude, which is our "input" Pomeron, also turns out to be 
small , a ~ 0.08 GeV~^ . Therefore, even at energy y/s = 18000 GeV^ we obtain for the radius of 
the constituent quark Rquark ~ 0.3 — 0.4 /m , that is still in two times smaller than the radius of 
the proton. This again supports the conclusion, that our interactions at LHC energy still far from 
the "black" disk limit. 

We also calculated the survival probability factor for different energies and it turns out to be 
similar to the values obtained in model proposed in Ref. ^8|. The parameters of the "input" 
Pomeron in our paper and in Ref. |28j are close, in spite of the fact that we did not involve 
any non-perturbative physics in our calculations. The intercept of the "input" Pomeron, which we 
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obtained , is A ^ 0.12 . This intercept is larger than the intercept of the paper l28|. The slope, 
a ~ 0.08 Gel^"^, is very close to the one of the paper j2Hl- Considering the impact parameter 
dependence of the amplitude of the exclusive central diffraction "hard" process, see Fig. UHl we see, 
that the maximum of the amplitude locates at peripheral impact parameter, 6 3 GeV~^ at energy 
of LHC The probability of this process at zero impact parameter is almost zero. 

The values of the parameters of our "input" Pomeron, the slope [pip ~ O.OSGeV"^ and the 
intercept A 0.2, give rise to the idea that this Pomeron is not "soft" (see Ref. [23j for the typical 
parameters of the soft Pomeron). The parameters that we obtained related to so called "hard" 
Pomeron. Therefore, the one of the result of this paper is the idea that, perhaps, we do not need to 
introduce "soft" Pomeron in order to describe the "soft" data. Considering the internal structure of 
the colliding hadrons as well as the unitarization corrections for the amplitude we can describe p-p 
data using only one, "hard" Pomeron. 
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Appendix A: 



In this appendix, as an example of our calculations, we calculate two diagrams, which contribute to 
the elastic amplitude: the first diagram of Fig. UHl and the first diagram of Fig. UHl 



1. 






4. 




Figure 15: Four diagrams of the three quark pair interaction. The wave line describes 
the eikonal amplitude for quark-quark interaction P{Y, b). 

We begin with the first diagram of Fig. ^1 There are six diagrams of this type and we need two 
vertices (see Eq. ^ and Eq. (jH}) for the calculations: 

and 

where qi is the momentum transferred of each single Pomeron, and k = qi + q2 + ■ So, we have 
for our diagram: 

D,,s,aak) = 6 / n ^n-,fe)e-^^/(^")e-^(^-^^3)^-^('-^'/^)^ {^tt' 6\k-q,-q,-q,)) (A.l) 

i=l ^ 

where Pq^g{qi) is the Fourier transform of amplitude given by Eq. Q. 

Since we have no simple analytic expression for Pg_g{qi) , we make Fourier transform for each 
function Pq_g{qi) and rewrite our expression in terms of Pg^gipi) : 

P,_,(g,) = j dH,Pg^g{k)e"^^'^. (A.2) 

Here bi is the impact parameter variable, which is conjugated to momentum qi of the single Pomeron. 
We also make Fourier transform from k to impact parameter variable b . We obtain: 

D,,spaUb) = 7^ / n d'kPg^gik) Jl[d'q, J dH e-''' (A.3) 

1=1 i=l 

Ell ^rb. g-fcV(6a) g-j^(,-,-,-3)^-^(,i-A;/3)^ _ g-. _ _ ) _ 

To make this expression easier for calculations we introduce new variables: 

2^1 = g2 - gs ; X2 = qi- k/3 ; X3 = qs , (A. 5) 

qi=X2 + k/3 ; q2 = Xi+x^; qs = . 
The Jacobian of this transformations is equal to unity; and we obtain: 
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-Dl^Spairs (^) 



6 



64 7r6 



n dH.P.^M / n d 



i=l 



1=1 

J_ 1,2 3_ ^2 _ _3_ ,.2 



d"^ k e* ''I (^2+fc/3) 



gif>2(^3+X-l)g*b3X-3 e-6^fc -8^^1-8^^'2 (2k / 3 - Xl - X2 - 2X3) . 

Performing the 0:3 integration we have: 



(A.6) 
(A.7) 



6 



4 647r6 



n d'b.P, 



i=l 



i=l 



gixi (62/2-63/2) giS2 (61-62/2-63/2) g-J-fc2-^x2-^x2 



Hc^'x, / (^2^ e^fcM+f'i/3+62/3+63/3) 

(A.9) 



Performing one integration over Xi and /c , we obtain the final answer for this diagram: 



J .^^ 



(A.IO) 



As the second example of the calculations technique we calculate the first diagram of Fig. IT?)1 
The main steps in calculation of this diagram are the same as in the calculation of the first diagram 
of Fig. ^1 Therefore, now we are not focused on the detailed explanations of main steps, but we 
rather clarify the principal points of the calculations technique. 

1. 







Figure 16: Four diagrams of the six quark pairs interactions. 
Two vertices of the Pomerons- quarks coupling have the following form (see Eq. ((7j) and Eq. 

Vup{k) 



e~6^ (5i+'?2+93)^-6^ (qi+gs+qe)^ {qi+q2+q3) (qi+gs+qe) 
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and 

Vdown{k) = e~8^(5i+'?4"'?2-95)^-6^(93+g6-gi/2-g2/2-g4/2-(?5/2)V, _ 

Here, the total momentum transferred is k = I]^=i Qi- We change the variables, in order to obtain 
a simple Gaussian expression for these vertices: 

xi - X2 = qi + q2 + qs Xi + X2 = q^ + q5 + qe ] (A.H) 
X3 = qi + qA-q2-q5 x^ = 3q-i/2 + 3qQ/2- k/2; (A.12) 
X5 = q5 xe = qe. (A. 13) 

And for old variables qi and k we have: 

k = 2xi gi = — xi/3 — X2 + X3/2 — X4/3 + X5 + ; (A. 14) 
q2 = 2xi/3-X3/2-x~A/3-X5 qs = 2x^/3 - Xq + 2xi/3 ; (A.15) 
q4 = Xi + X2-X5-X6 q^ = xr,qQ = x^. (A.16) 

Here momenta qi are not independent. For these momenta we used the delta function constraint: 

qi = k - q2 - qz - qi - q^ - q&- 

Using Eq. ()A.14|1 we calculate the Jacobian of the variable change, it is equal to | . In new variables 
the vertices look very simple: 



up 

and 



Kp(fc) = e-^"''-^"', (A.17) 



1 ^2 1 ^2 



^do«,n(A;)=e-^^3-67j-4. (A.18) 

We consider the exponents that stem from the Fourier transform from momentum to impact 
parameter representation: 

e-^^^e^nli9^»^^« , (A. 19) 

where we also have 

n p,-M ^ n ^.-.(^0 ■ (A.20) 

2=1 i=l 

Putting in Eq. ()A.19|) substitutions of Eq. ()A.14|) we obtain for the exponents of Eq. ()A.19|) : 

gixi{-26+262/3+2t;3/3+fe4-fci/3) gjx2(fe4-fei) gi X3(bl/2-b2/2) gix4(-b2/3+2fe3/3-bl/3) (^p^ 21) 

gix5(-62-f'4+f'5+f'i) Xfi{-b-j.-b4+br,+b-i) 22) 

The vertices of Eq. ()A.17|) and Eq. ()A.18|) have no dependence on x^ and Xg variables. Therefore, 
in integration over X5 and x^ we obtain the following delta functions: 

+ 61 - 62 - ^4) ^ h = b2 + h-bi , (A.23) 

and 

+ 61 - &3 - &4) ^ K = h + b,-h. (A.24) 
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In t his ca se, after the integration over 65 and , the answer for the functions Pq-q{hi) ( see 
Eq. (IA.20jl 'l: looks as follows 

6 4 

n P.-Mi) ^ n Pi-<i(bi) ^.-.(^2 + h- h) Pq-q{h + K- hi) , (A.25) 

i=l i=l 

Of c ourse, if we need, w e should also replace 65 and in the other exponents of Eq. (jA.21|l substituting 
Eq. ()A.23j) -Eq. ()A.24|l .In our particular diagram this replacement has no place since exponents of 
Eq. ()A.21|1 do not depend on 65 and b^. 

We are ready to write the answer for our diagram. Performing a simple integration over variables 
Xi , where i = 1 — 4 , we ob tain (we use Gaussian functions of Eq. ()A.17|) and Eq. ()A.18|) with the 
four exponents of Eq. ()A.21|) ): 



Di,ePom{b) = l[d^b^ Pcj-M P,-,{b2 + &4 - bi) Pq^qibs + 64 - &i) (A.26) 

^ i=l 

(26-2b2/3-2fe3/3-64+6i/3)2 - f (h-hf - f (fei-fea)' - ^ (63-&2/2-bi/2)2 

All other diagrams, which contribute to the elastic amplitude, we calculate using the same meth- 
ods which have been described above. The answer for entire amplitude, which included expression 
for all possible diagrams of our model is written in Appendix B. 

Appendix B: 

In this appendix we present the resulting expressions for the full elastic amplitude, order by order. 
We have for the amplitude ( see Eq. 0): 

A{S, b) = Aipair{s, b) - A2paivs{s, b) + A'ipairsis, 6) + ... + Agpairs{s, b) . (B.l) 

The answer for the contribution to the amplitude from the one pair of quarks , Aipair{s., b) , is given 
in Eq. (fTT|) . So, we start from the two pairs contribution. The diagrams of this contribution are 
shown in Fig. El 



1. 2. 




Figure 17: Two diagrams for the interactions of two quark pairs. 
We have for A{s, b)2pair '■ 

A^paUsM = ^ /nrf^6.P,-,(6.)e^^(^^/^+^V^-^>-t(^-^^)^ + (B.2) 

97r„2 /. 2 ^ ^ ^ ^ ^ 

+ ^ JUd'h P,-M e-'-^'^'^^'^l'-'^'' - f (^-^^)^ (B.3) 

^ 4 = 1 
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The diagrams for the three pairs of interacting quarks are shown in Fig. ^1 We have for this 
contribution: 



^ / f[dH,P,.,{b,)e-'^^'-'^/'-'-^^^^^^^ (B.4) 



+ 
+ 
+ 



TT 

9 a 



(B.5) 



i=l 
3 „ 2 



n h Pg-M P<l-,i^b -h- &l)e-t (3^-3^-1/2-3^/2)2 - f ^g.G) 



i=l 



27a3 - 3 



-3a (b-fe2/2-63/2)'^ - f (b2-bi)" - f (bi-bs) 



For four pairs amphtude we have the diagrams of Fig. ^1 

1. 2. 



3. 




Figure 18: Five diagrams for interactions of four quark pairs. 
For this contribution we have: 



(B.7) 



+ 



IT" 



i=l 



-f(6fe+6i-2fe2-263-3fe4)2-f(6i-b4)2-f (61-62)2-^(61/3+62/3-263/3)2 



+ 



i=l 



g-f (66-61-262-63-264)2-1(61-63)2-1(61-62)2-41^(61/2-62/2-63+64)2 _^ 



4 ^ 3 



+ 1^ / n P,-,ibi -b2-h 

g-^(26-62-63)2 -f(6i-63)2- f (61-62)2 ^ 



(B.8) 

(B.9) 
(B.IO) 

(B.ll) 
(B.12) 

(B.13) 
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e 4 
9a 



i=l 



}{bl-b2? 



+ 



i=l 

^-^{4b-bi-b2-h-b4)^ - f (64-63)2 - f (61-62)2 - i|(6i+62 63-64)2 



The five pairs interaction diagrams are shown in Fig. ^1 

1. 2. 3. 




Figure 19: Five diagrams for the interactions of five quark pairs. 
In this case we have: 



(B.14) 

(B.15) 
(B.16) 

(B.17) 



A 



bpain 



[S,b) 



r ^ 

i=l 



e 2 
9a4 



a^^(2b+6i/3-2b2/3-2fe3/3-64)2-f (61-64)2- 1(61-62)2-^(61+62-63)2 _^ 



+ -;^JIld'b^ Pg-M P,_,(36 - 61 + 62 - &3 - &4) 



i=l 

g-^(36/2-6i/2+fe2-363/2-64/2)2-f (62-64)2-1(36-61-63-64)2-^(61-62/2-64/2)2 



+ ^ lid'h P,-M P,-,(36 - 62 - h) 
^ i=i 



2^^(36/2-61/2-62-63/2-64)2-1(61-64)2-1(61-63)2-^(61/2+62-63-64/2)2 _^ 



+ -^JUd'h P,-M P,_,(36 - 62 - 64) 



i=l 

g- ^ (36-61-62/2-63/2- 



)2-f (61-62)2- 1(62-63)2-^(61/2-62/2+63-64)2 _j_ 



+ / n rf' fe. P,-M Pg-,{Sb - 62 - &3) - &2 - &3 

i=l 



(B.I8) 

(B.19) 
(B.20) 

(B.21) 
(B.22) 

(B.23) 
(B.24) 

(B.25) 
(B.26) 
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^_3^(26-fe2-fe3)2-f (?i-fe2)2-f(?2-?3)2-f(6l-fe/2-63/2)2 



(B.27) 



The diagrams with the interaction of six quark pairs are shown in Fig. ^1 For this contribution 
we have: 



A, 



Gpairs 



6 



TT' 



j \{<fh P,~,{h) P,-,(&2 + &4 - h) P,^,{h + h-h 

i=l 



^(2fe-262/3-2b3/3-64+fei/3)2 - 1(64-61)2 - f (61-62)2 - 2ji(53_52/2-6i/2)2 _^ 



in r 

+ J lid'b, P,_^{bi) Pg^g{3b -h- b,) P,_g{b2 + 64 - h) 

^ i=l 

(86-61/2-62+63-364/2)2 - f (62-61)2 - f (61-64)2 - ^ (63-61-62)2 _^ 

36 a'^ /• 4 

+ ^ / n ^5-<?(3& + &1 - ^2 - &3 - &4) n-g(&2 + &4 - &1 

^ i=l 

(36+61/2-62-63-364/2)2 - f (62-61)2 - f (61-64)2 - ^ (63-61/2-62/2)2 _^ 
c 4 4 

+ ^/ n ^'^'-P<?-5(&i)n-g(3&-^2-&4)P,-,(3&-6i-63) 

^ 4 = 1 

^-^ (36-61-62/2-63/2-64)2 - f (62-61)2 - f (62-63)2 - 2ji (64-63-61/2+62/2)2 

The diagrams with the interaction of seven quark pairs are shown in Fig. EOl 

1. 2. 
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(B.28) 

(B.29) 
(B.30) 

(B.31) 
)(B.32) 

(B.33) 
(B.34) 

(B.35) 



Figure 20: Two diagrams for the interactions of seven quark pairs. 
For these diagrams we have : 



A 



Ipairs 



IT" 



i=l 



(B.36) 



+ 



P,_,(26 - 64 - b;) P,_,(62 + &4 - bi) P,_,(6i - 62 + &3) (B.37) 

^-^(36-61/2-62/2-63-64)2 - f (64-61)2 - f (61-62)2 - 2f{b:i~^b2+bi/2-b4/2f (^g gg^i 
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n dH,Pg.,{k) 



i=l 



P,-,(36 - 64 - &3) P,-,i3b - 64 - &2) P. 



«-9l"2 



61+64 



g-^(36+36i/2-62-63-364/2)2 - f (64-61)2 - f (61-62)2 - ^(b-i-b2/2-bi/2)^ 



(B.39) 

(B.40) 
(B.41) 
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8 interacting pairs 
of quarks 

It ' 



9 interacting pairs 
of quarks 




Figure 21: Two diagrams for the interactions of eight and nine quark pairs. 

Finally, there are only one type of diagrams with the eight and nine quark pairs interactions, see 
Fig. EH 

So we have for the diagrams with the interaction of eight quark pairs: 



a 



J l[dH, P,.,{h) P,_,(26 - h) P,_,(26 - 62) 



i=l 



63) Pg-g{bi -62 + ^4 



-^ib-bi/2+b2/2-b4f - f (fe2-fei)2 - f (2fe-6i-62)2 - ^(26-61/2-^62/2-63)' 



and the diagram for nine quark pair interactions gives: 



^9pairs{Sj b) 



a 

9^ 



r ^ 

J l[dH, P,_,(6,) P,_,(26 - 61) P,_,(26 - 64) 

i=l 



Pg-,i2b - 61 



bi)Pq^q{b^ - 62 + bA)Pq-q{bl - 62 - &3 



g-^(6+6i/2-63-64/2)2 - f (64-61)2 - 1(26-61-64)2 - ^(61/2-62-64/2)2 



(B.42) 

(B.43) 
(B.44) 



(B.45) 

(B.46) 
(B.47) 



In our numerical calculations we used only the A{s,b)ipair-A{s,b)j pairs terms of the amplitude. 
The term A{s, &)6pajrs-^(s, b)-jpaiTs works only at energies close to the LHC one. 



Appendix C: 

As example of expression for the calculation of the survival probability (SP) in the framework of our 
approach we consider the term of elastic amplitude with three interacting quark pairs: 



A{sM,,a^rs = ^^J\{dH.-,P,,,{h:)e'''f^^^^^^^^^ 
^ i=l 

+ '^J^ jj^Sh P,-M e-t (3^-^-^2-53)' - f (^2-6.)' - ¥ (b-bsf + (c_2) 

^ i=l 

+ % fUd'h P,-MP,-,iSb - 62 - 6,)e-f (3^^"36i/2-362/2)2-f (6,-6i)2 ^ ^^ 3) 
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3 „ 3 



^ i=i 

Let us rewrite the first term of tliis expression. Using receipt of the Eq. (jSSI), we rewrite: 



l[P<^-M - F'sp^3,arrsibl,b2,h) = P,^,{h) P,^M P,^M + (C.5) 



The expression of the r.h.s. with Eq. ()(7.5|1 instead of IlLi Pq-qipi) in the first term of elastic 
amphtude. Eq. ()C.1|) gives the answer for the first term of SP amphtude A2,pairs{.s,h). Of course, 
obtaining A(s, h)^pairs we must to perform such replacement in each term of Eq. ()C.1|) . Doing so we 
obtain: 

A{s,h)^pairs = (C.7) 

6a^ " ^ 



27r3 



i=l 

+ ^ / n rf^&.i^|p,3.a..(&l,&2,63)e-f («-'^-'^3)-f + 
i=l 

+ ^ / n Flp,s,aabi, b,) PU^b - 6, - W)e-f (3^-3V2-3?./2)^ - f (^.-M^ + 

i=l 

+ ^ / n ^^^'.i^.W.(&i'^2'&3)e-^"(^-^^/-^^3/2)^-^ . 
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